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Recent evidence for direct CP violation in non-leptonic charm decays cannot be easily accom- 
modated within the Standard Model (SM). On the other hand, it fits well in new physics models 
generating CP violating AC = 1 chromomagnetic dipole operators. We show that in these frame- 
works sizable direct CP asymmetries in radiative D — s> P + P~j decays (P = it, K), with Mpp close 
to the p or the (j> peak, can be expected. Enhanced matrix elements of the electromagnetic dipole 
operators can partly compensate the long distance dominance in these decays, leading to CP asym- 
metries of the order of several percent. If observed at this level, these would provide a clean signal 
of physics beyond the SM and of new dynamics associated to dipole operators. We briefly comment 
on related CP violating observables accessible via time dependent D(D) —¥ P + P~"/ studies and 
angular decay product distributions in rare semileptonic D decays. 
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I. INTRODUCTION 

A significant evidence for direct CP violation in D — > 
P + P~ decays (P = ir, K) has recently been reported 
by the LHCb [1] and by the CDF [2] collaborations. 
Both experiments find a non- vanishing value for Aac p = 

o-k+k- — a-TT+TT-i where 



_ T(D° -> /) - r(L>° -> /) 
T{D° ->■ /) +r(D° ->■ /) 



(1) 



Combining these recent results with older measurements 
leads to the following world average [2] 



Acxp 
a CP 



-(0.67 ±0.16)% 



that differs from zero by about 4<r. 

The theoretical interpretation of this result is puz- 
zling. The value in Eq. (2) exceeds by a factor 5-10 what 
is naturally expected in the Standard Model (SM) (see 
e.g. Rcf. [3], and the more recent analyses in Ref. [4, 5]). 
However, wc cannot exclude that such result has a SM 
explanation due to the non-perturbativc enhancement of 
penguin- type hadronic matrix elements [6-8]. On the 
other hand, this value can naturally be accommodated 
in well-motivated extensions of the SM. In particular, it 
fits well in models generating at short distances a sizable 
CP violating phase for the effective AC = 1 chromomag- 
netic operators [3, 4, 9, 10]. 

Given this situation, it is important to identify possi- 
ble future experimental tests able to distinguish standard 
vs. non-standard explanations of Aacp- An interesting 
strategy that makes use of CP asymmetries in various 
hadronic D decays (necessarily including neutral mesons) 
has recently been proposed in Ref. [11]. However, this 
strategy is effective in isolating possible non-standard 
contributions to Aacp only if they are generated by ef- 
fective operators with a A J = 3/2 isospin structure. This 
is not the case for the well-motivated scenario with a new 



CP violating phase in the AC = 1 chromomagnetic oper- 
ator. As we point out here, in the latter case an efficient 
strategy is obtained by measuring CP asymmetries in ra- 
diative D decays. 



II. SHORT-DISTANCE EFFECTIVE 
HAMILTONIAN 

The first key ingredient of our strategy is the strong 
link between the AC = 1 chromomagnetic operator, 



u L a^T a g s G^c R , 



(2) and the AC = 1 clcctromagnetic-dipole operator, 



(3) 



(4) 



In most explicit new-physics models the short-distance 
Wilson coefficients of these two operators (C73) are ex- 
pected to be similar. Moreover, even assuming that only 
a non- vanishing C 8 is generated at some high scale, the 
mixing of the two operators under the QCD rcnormaliza- 
tion group (RG) implies C7 8 of comparable size at the 
charm scale. The same is true for the pair of operators 
with opposite chirality Q' 7 g , obtained from with the 
replacement L -h- R. 

To quantify the size of these coefficients, we normalize 
the effective Hamiltonian describing the AC = 1 new- 
physics contributions as 



H 



cff-NP 
|Ac|=l 



G 



V2 



(5) 



The complete list of potentially relevant operators can 
be found in Ref. [4]; however, for the purpose of our 
analysis we can restrict our attention only to Q.7,8 and 
<2' 7 g . Assuming the initial conditions of these operators 
are generated at some scale M > m t , taking into account 



2 



the RG evolution of the operators at the leading log level 
(assuming only SM degrees of freedom below the scale 
M), leads to [12] 



Cf (m c ) = V vdf'iM) + 8 (t? - 1) C ( 8 ''(M) 



0/ 



C ( 8 '\m c ) 



where 



and 



n 



~a s (M)~ 


21 


~a s (m t y 


23 


"a s (?7i fc )" 


_a s (m t ) 




_a s (m b )_ 




_a s (m c ) 


~a s {M)~ 
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21 


~a s {m t )~ 
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a s (m b ) 


_a s (m t )_ 




_a s (m b )_ 




_a s (m c ) 



(6) 
(7) 

(8) 
(9) 



Following the analysis in Rcf . [9] , the new-physics con- 
tribution to Aacp induced by Qg, can be written as 



|Aag£| 



-1.8|Im[Cr(m c )]| , 



(10) 



where the numerical value assumes maximal strong 
phases and is affected by 0(1) uncertainties due the the- 
oretical error on (PP\Q$\D). Assuming this contribu- 
tion saturates the experimental value of Aacp leads to 
|Im[C| IP (TO c )]| ~ 0.4 x 10 -2 . If we further assume that 
the initial scale M is around 1 TeV, and that at this scale 
|Cf p (M)| < \Cf p (M)\, the RG evolution implies 



|Im[Cr( 



|Im[C 8 NF (m c )]| «0.4x 10" 



(11) 



This is for instance what happens in supcrsymmetry, 
where the gluino-mcdiated amplitude proportional to 
(5® R )i2 leads to the initial condition 



£>SUSY 



(msusv) = (4/15)C 8 SUSY (m SUSY ) . (12) 



Taking into account the 0(1) uncertainties in the deter- 
mination of |Im[C 8 N[p (m c )] | , and the additional uncertain- 
ties in the initial conditions of Cf p (M), we consider the 
following range for lm(Cf F ) at the charm scale 



|Im[C? p (m c )]| = (0.2 - 0.8) x 10" 



(13) 



The same range holds for Im(C' 7 ), if the leading contri- 
bution to Aacp is generated by Q' s rather than Qg . 

At low energies CV receives contributions also from the 
mixing with the SM four-fermion operators. However, to 
a good accuracy these contributions are CP conserving. 
The leading effect is the two-loop mixing between C7 
and Ci'2 [13]. According to the analysis in Ref. [13], 
integrating out also light quark loops one obtains 



SM-cff 



(m c )| = (0.5 ±0.1) x 10" 



(14) 



with an 0(1) strong phase and a negligible CP-violating 
phase (more than two orders of magnitude smaller) . 

If the contributions in Eqs. (11) and (14) where the 
dominant contributions to radiative D decays, we could 
expect 0(1) direct CP asymmetries in these modes. As 
we discuss below, this is not the case due to genuine long- 
distance contributions that dominate the decay rates. 



III. SHORT- VS. LONG-DISTANCE 
CONTRIBUTIONS IN D — > Vj 

The second important ingredient of our analysis is the 
observation that in the Cabibbo-suppresscd D —> U7 
decays, where V is a light vector meson with uu va- 
lence quarks (V = p°,uj), Q7 and Q' 7 have a sizable 
hadronic matrix element. More explicitly, the short- 
distance contribution induced by Q^p , relative to the to- 
tal (long-distance) amplitude, is substantially larger with 

respect to the corresponding relative weight of in 
D -> P + P~ decays. 

The decay amplitudes for D —> V7 decays can be de- 
composed as follows 



A[D(p)^V(p,ih(q,e)] 



-A^ v [(e*q)(e*p)-(qp)(te*)} 



(15) 



with the corresponding rates 



32tt 



2 \ 3 



[\A PV \ 2 + \A PC \ 2 ] 



(16) 

The short-distance contribution induced by Q7 to the 
effective couplings Ap V PC is 



(^pc(PV)) s ' d ' 



eQ u G F rn c 
v/2 2tt 



CV(m c 



J l(2)' 



(17) 



where TXp) are defined by 



{V{p,e)\uq v a^(l + 75 )c\D(p}) = -2ie^r a ^p"TY 
+T 2 v [(m 2 D -m 2 v )i*^(rp)(p + pr] , (18) 



and TY 



rpV 

1 2 



rpV 

1 {D) 



a/9- 



T, 



via the identity 75 
A recent sum-rule estimate finds [14] 



(D) 



0.70(7). We note in passing that at leading or- 
der in a s and in the infinite charm quark mass limit, 
heavy quark symmetry predicts TL = V^(0)/(1 — 

mv/mu) [15], where (q 2 = (p — p) 2 ) 

V v (q 2 ) 

(V(p, e)\u^c\D{p)) = 2i {D) " e^ af3 p„p a ~e p . (19) 
mo + my 

This matrix element enters semileptonic D — > V decays 
and thus ^(0) can be accessed experimentally. Un- 
fortunately for the interesting D — > (p,uj)£v transitions, 
no such analyses are available at present. On the other 
hand, in the heavy charm quark limit, T^ can be related 
to hadronic matrix elements entering radiative B — > V 
transitions - TY B y Starting from the (quenched) Lattice 
QCD estimate T ( P B) = 0.20(4) [16], and running it in both 
the perturbative matching scale (from /i b = 4.6 GeV to 
/x c = 1.4 GeV) as well as the heavy quark mass scaling 
(including leading power corrections) [16, 17], we obtain 
T? D , w 0.7(2). Using instead existing sum-rule estimates 



3 



of [18] typically leads to 0(20%) larger values. Con- 
sequently we employ the value of TY D -, with a conservative 
uncertainty estimate of 



rpp 

1 (D) 



t; 



(D) 



0.7 ±0.2, 



which leads to 



( AP>u \s.d. 
^PC,PVJ 



0.6(2) x 10- 



C 7 (m c ) 



0.4 x 10" 2 



(20) 



(21) 



The contribution induced by Q' 7 is obtained with the re- 
placement CV — > ±C 7 in -^p'c(pv)- 

The only D° V°j decays observed so far are the 
K* and <\> modes [19]. The observed rates satisfy to 
a good accuracy the relation B(D — s- K*j)/B(D —> 
K*p°) = B(D <fo)/B{D -> (jyp ), generally ex- 
pected by vector meson dominance (VMD). The three 
Cabbibo-suppressed D° — > V°j modes, V° = p° 7 uj 7 (f>, 
are expected to have similar rates. 1 We can thus esti- 
mate the typical size of their long-distance amplitudes 
\(A^ c y- d -\ ~ \(A% y-*-\ as follows 



K^pc,pv) Ld ' 




32tt 








m D 





-3 



T(D -> Vj) 



5.8(4) x 10" 



for V 



1/2 



(22) 



In the limit where the strong phases of the amplitudes 
have a mild energy dependence, and assuming we can 
neglect the weak phase of the long-distance amplitude 
(see Sec. V), the direct CP asymmetry, defined in Eq. (1), 
can be decomposed as 



where 



\ a v~f \ = 2 Cwoak I sin(A^ strong )| , 
|lm(^ CiPV ) 8 ' d -| 



C 



weak 



(^pc : pv) 1,d ' 



(23) 



(24) 



As a result, according to Eqs. (21) and (22), in the p and 
ui modes the CP violating asymmetries can reach 10% 
for maximal strong phases: 



-- 0.04(1) 

10~ 5 



Im[C 7 (m c )] 



0.4 x lO- 2 

1/2 

< io% 



(25) 



1 According to explicit VMD predictions [20] B(D — > fry) and 
B(D — > wy) are very similar, possibly a factor ~ 2 smaller than 
B(D° -> <jry). In the following we assume 3{D —¥ (p, uj)^y) > 
1CT 5 . 



The case of the <f> resonance, or better the \K + K~ r y) fi- 
nal state with Mkk close to the 4> peak, is more involved 
since the hadronic matrix element (18) vanishes, in the 
large m c limit, if V is a pure ss state. However, as we dis- 
cuss in more detail in the next section, a non-negligible 
CP asymmetry can be expected also in this case for two 
main reasons: 1) the matrix element in (18) is not iden- 
tically zero even for V = <f), both because 0(Aqcd /m c ) 
corrections and because of the tiny uu component of <fi; 
2) non-resonant contributions due to (off-shell) p and ui 
exchange can also contribute to the \K + K~^) final state. 



IV. THE D -> K+K-j CASE 

The decay amplitudes for D — > P + P~j decays can be 
decomposed in full generality as follows 



A[D( P ) -> P H 

-iM(s,u) e^ Q(9 <A*V 
+E(s,v) e*jq»(qp+-qp-) 



(p+)p-(p_)7(g,e)] = 
P a (p+-P-Y j 

qp(p+ 



P-) 



(26) 



where s = (p+ + p_) 2 and v = {qp+ — qp-)- In the limit 
where we consider at most electric and magnetic dipole 
transitions (or neglecting higher order multipoles), we 
can neglect the v dependence of the form factors. In 
this approximation, the differential rate as a function of 
s 



Mp P can be written as 



ds 



32tt 



D 



[\M(s)\ 2 + \E(s)f 



(27) 

where T (s) = y/s(l - Amp/ s) 3 / 2 / (48tt). 

If the amplitude is dominated by the exchange of vector 
resonances we can decompose M and E as follows 



M(s) = 



9 V pp 



^PC 



E(s)=J2 



9pp 



Ml 



(28) 
(29) 



where gp P is the V —> PP coupling, defined such that 
r(V PP) = gppT (Ml). It is then easy to check that 
in the limit of a single narrow resonance, integrating over 
s, we recover T(D -> PPj) = T(D -> Vj)xB(V -> PP). 

In order to estimate the maximal direct CP asymmetry 
in the D — > K + K~~y case, with Mkk close to the </> 
peak, we evaluate M(s) and E(s) summing over the three 
light vector resonances (V = p,ui,(f>) with the following 
assumptions: 

• In all cases we use the parametric form in Eq. (22) 
to estimate the overall magnitude of ^4pc(pv)' as " 
suming further B(D — > (p,uj)j) > 10~ 5 . 

• For V = p, u> we assume the weak phase of ^4pc(pv) 
is Cweak, while for V — 4> we use rC, W cak- Here r — 



4 



0.3(1) is the typical annihilation suppression factor 
in non-leptonic D decay amplitudes [7, 21], that we 
apply to the the matrix element in Eq. (18) in the 
V = 4> case. 

• For V — p,u> we the effective coupling to K + K~ 
to g]^+K- — 3, as expected by SU(3) symmetry 
given that gP^ ~ 6. 

Under these hypotheses, and assuming maximal and 
smoothly varying strong phases for the contributions 
with different weak phases, we find 

|aK+A:- 7 | max « 2% , 2m K < yG < 1.05 GeV , 
\a K +K-y\ max « 6% , 1.05 GeV < *J~s < 1.20 GeV . 

(30) 

In the first bin, close to the cf> peak, the leading contri- 
bution is due to the ^-exchange amplitude. The con- 
tribution due to the non-resonant amplitudes plays a 
significant role far enough from the <f> peak, where the 
charge asymmetry can become larger. However, it must 
be stressed that away from the </> peak the overall rate of 
the D — > K + K~^f process is significantly reduced. 

V. DISCUSSION 

In order to establish the significance of these results, 
two important issues have to be clarified: 1) the size of 
the CP asymmetries within the SM, 2) the role of the 
strong phases. 

As far as the SM contribution is concerned, we first no- 
tice that short-distance contributions generated by the 
operator Q7 are safely negligible: using the result in 
Ref. [13] we find asymmetries below the 0.1% level. The 
dominant SM contribution is expected from the lead- 
ing non-leptonic four-quark operators, for which we can 
apply the general arguments presented in [4]. The CP 
asymmetries can be decomposed as 

|of*| « 2£ lm(R s f M ) « 0.13% x lm(R s f M ) , (31) 

where £ = iV^V^hl/IVcsV^sl and R^ M is a ratio of sup- 
pressed over leading hadronic amplitudes, naturally ex- 
pected to be smaller than 1. This decomposition holds 
both for the / = ttit, KK channels discussed in Ref. [4] 
and for the / = U7 case analyzed here. The SM model 
explanations of the result in Eq. (2) require R^ KK ~ 3. 
While we cannot exclude this possibility from first princi- 
ples, a further enhancement of one order of magnitude in 
the D —> V7 mode is beyond any reasonable explanation 
in QCD. As a result, an observation of |ay 7 | > 3% would 



be a clear signal of physics beyond the SM, and a clean 
indication of new CP-violating dynamics associated to 
dipole operators. 

Having clarified that large values of |ay 7 | would be a 
clear footprint of non-standard dipole operators, we can 
ask the question if potential tight limits on |ay 7 | could 
exclude this non-standard framework. Unfortunately, 
the uncertainty on the strong phases does not allow 
to draw this conclusion. We recall that the maximal 
values in Eqs. (25) and (30) can be reached only in 
the limit of maximal constructive interference (namely 
of ±7r/2 strong phase difference) of the amplitudes 
with different weak phases. The calculation of light- 
quark loop contributions in Ref. [13] does suggest the 
presence of large strong phases in these amplitudes; 
however, we cannot exclude destructive interference 
effects leading to \ayy\ = O(0.1%) even in presence 
of a non-standard CP-violating phase in the dipole 
operator. In principle, this problem could be overcome 
via time-dependent studies of D(D) Vj decays or 
using photon polarization, accessible via lepton pair 
conversion in D —> V(j* —> £ + £~); however, these types 
of measurements are certainly more challenging from the 
experimental point of view. 



VI. CONCLUSIONS 

Radiative D — > P + P~j decays, with Mpp close to 
the p or the <j> peak (for P = it or K, respectively), could 
help to shed light on the origin of CP violation in the 
charm system. If the experimental result in Eq. (2) is due 
to non-standard dynamics involving dipole operators, we 
can expect significantly larger direct CP asymmetries in 
these radiative modes. As we have shown, evidence of 
|app 7 | > 3% would be a clear signal of physics beyond the 
SM, and a clean indication of new CP-violating dynamics 
associated to dipole operators. 
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